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Nanoscale and microscale confinement of biopolymers naturally occurs in cells and has been re- 
cently achieved in artificial structures designed for nanotechnological applications. Here, we present 
an extensive theoretical investigation of the conformations and shape of a biopolymer with varying 
stiffness confined to a narrow channel. Combining scaling arguments, analytical calculations, and 
Monte Carlo simulations, we identify various scaling regimes where master curves quantify the func- 
tional dependence of the polymer conformations on the chain stiffness and strength of confinement. 
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What is the effect of confinement on the shape of a 
biopolymer? With recent advances in visualizing and 
manipulating macromolecules on ever shrinking length 
scales, an answer to this question has gained increasing 
importance. In the crowded environment of a cell the 
conformations of cytoskclctal filaments are highly con- 
strained by other neighboring macromolecules. This con- 
finement largely alters the viscoelastic response of entan- 
gled biopolymer solutions There is growing interest 
in manufacturing nanostructures such as nanopores Q 
and nanochannels 0, Q for investigating and manipu- 
lating DNA with improved technologies aiming towards 
smaller and smaller structures. Hence an improved un- 
derstanding of the effect of confinement on biopolymer 
conformations has potential implications for the design of 
nanoscale devices in biotcchnological applications. Simi- 
larly, microfluidic devices have been used to explore con- 
finement effects on actin filaments and DNA [1,0] ■ What 
makes the confinement of biopolymers both a challenging 
and interesting problem is that biopolymers, unlike their 
synthetic counterparts, arc generally stiff on a length 
scale much larger than their monomer size. The persis- 
tence length t p , the scale below which bending energy 
dominates over thermal fluctuations, is approximately 
50 nm for DNA [| and 16 ^m for F-actin [9j. Depending 
on whether the contour length L is smaller or larger than 
the persistence length we may distinguish between stiff 
and flexible chains. 

For cellular systems as well as for nanoscale devices, 
biopolymers are confined on length scales comparable 
with their persistence length £ p such that the polymer's 
intrinsic bending stiffness plays a decisive role for its con- 
formations. For simplicity, consider a cylindrical tube 
of diameter d. Upon balancing the bending stiffness of 
a chain with thermal energy, Odijk [Io| has identified a 
length Ld measuring the typical distance between succes- 



sive deflections of the chain, L\ 



d 2 £ p ; see Fig.ffl This 



suggests to use the number of collisions c = L/L^ per 
filament length L as a natural dimensionless parameter 
to measure the strength of confinement and e = L/£ p to 
measure the flexibility of a polymer. The physics in the 



strong confinement regime (c> 1) is genuinely different 
from the regime where the radius of gyration Rq of a 
long flexible chain (with L ^ £ p ) becomes comparable 
to d. In the latter case of weak confinement of a flexible 
chain the shape of the polymer is distorted due to self- 
avoidance between distant segments along the polymer 
chain. Then, according to de Gennes' blob picture [ill ], 
we may represent the conformation of the polymer as a 
linear chain of nonpenetrating spheres of radius d, where 
each sphere is described by Flory's theory. This picture 
results in the following scaling law for an extension of the 
polymer along the tube axis: 



i?H ~ L (£ p /d) 2 ^ . 



(1) 



For strong confinement the Odijk length becomes the 
analogue of the blob size, below which the polymer may 
be considered as free. Thus the fraction of contour length 
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FIG. 1: (Color online) Scaling regimes for confined biopoly- 
mer conformations as a function of polymer flexibility e = 
L/lp and confinement strength c = L/La- In the flexible 
regime (e 3> 1) two scaling regimes are known [dark grey 
(red)]: free coil behavior for weak confinement and de Gennes 
scaling for intermediate confinement. In the parameter range 
most relevant for biopolymers [light gray (green], one has to 
distinguish between weak confinement of stiff polymers and 
strong confinement for chains of arbitrary stiffness. 
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stored in thermal undulations decreases with decreasing 
tube dimensions [Tpj : 



(L-R\\)/L~L d /tp. 



(2) 



This scaling law should apply equally well for stiff and 
flexible chains as long as the collision parameter c is suf- 
ficiently large. For flexible chains, e> 1, this is the case 
if the deflection length Ld is less than the persistence 
length t p . For stiff chains, e <C 1, the deflection length 
has to become smaller than the total filament length be- 
fore there is any stretching. There is an additional regime 
of weak confinement (d < L < Ld) where the average ori- 
entation of the filament becomes aligned with the tube 
axis. These various scaling scenarios are summarized in 

Fig. hi 

The purpose of this Rapid Communication is to go be- 
yond this qualitative scaling picture and provide a quan- 
titative study of the conformations of biopolymers in con- 
fined geometry. Our focus is on the parameter range 
that is most relevant for cellular systems and nanoscale 
devices [light gray (green) region in Fig. [T] , where self- 
avoidance effects may safely be neglected. For specificity, 
we consider a wormlikc chain in a soft harmonic poten- 
tial of cylindrical symmetry and strength 7. Indeed, the 
harmonic potential is the simplest model to represent a 
tubclikc confinement and it has the advantage of being 
amenable to analytic calculations. Thus the Hamilto- 
nian for the contour r(s) parametrized in terms of the 
arc length s reads 



n[r(s)} = - 



ds 



d 2 r(s) 
ds 2 



dsvl(s), (3) 



where K = ipk^T is the bending stiffness, and = (x, y) 
are the components of the contour perpendicular to the 
tube axis. 

A variety of analytical results have been obtained so 
far for unconfined wormlike chains. For instance, the 
tangent-tangent correlation functions [l2T | and moments 
of the end-to-end distributions have been calculated ex- 
actly [ill, [HI . Further results like the probability distri- 
bution function of the end-to-end distance R have been 
calculated for stiff chains [l5| within the weakly bend- 
ing rod (WBR) approximation, where one considers only 
small transverse bending fluctuations with respect to a 
straight contour r(s) « (r±(s),s). For confined chains 
the WBR limit amounts to assume that the filament is 
aligned with the tube axis, and may be employed to study 
the asymptotic regime of strong confinement, Ld <C L. 
Then, using the equipartition theorem one obtains for 
the correlation function of the transverse undulations in 
Fourier space: 



(x(k)x(-k)) = 



c(fc 4 + 4L" 



(4) 



where Ld = (4K/7) 1 / 4 . This simple result forms the basis 
for all subsequent analytical calculations. For example, 
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FIG. 2: (Color online) Left: marginal probability distribu- 
tion function of polymer configurations in a plane containing 
the tube axis for c = 4 and e — 2. The probability den- 
sity increases from dark to light colors. Right: normalized 
tangent-tangent correlation function (j>( s ) = (t(s) ■ t(0.5)) for 
e — 0.1 and various confinements ranging from c = 3 to 
c = 20. MC simulations are represented by symbols, whereas 
the solid lines show the analytical approximations of <j)(s) as 
from Eq. (0 ■ Dashed line: exponential decay of an unconfined 
chain. 



it implies that the local transverse mean-square displace- 
ment is given by (r^) = L d /4£ p . 

In general, however, analytical calculations are not fea- 
sible. In order to investigate the full range of parameters 
we employed a standard Monte Carlo (MC) scheme us- 
ing a discretized version of the wormlike chain model. 
This allows us to go beyond the WBR limit and calcu- 
late a range of obscrvablcs which arc directly accessible 
in single-molecule experiments. The polymer was repre- 
sented by a chain of N segments tj approximating the 
continuous contour. The inextensibility constraint was 
imposed along the whole polymer by fixing the segments 
length to the value L /N. The cylindrical symmetric har- 
monic potential was calculated at the end points of each 
segment. 

During the simulation, both ends of the polymer were 
assumed to be completely free, in both position and 
orientation. The initial configuration was chosen suffi- 
ciently close to full stretching, in order to ensure a fast 
convergence of the MC algorithm. A new configuration 
was generated by changing the orientation of a randomly 
chosen segment and accepted according to the standard 
Metropolis algorithm. We have not considered effects re- 
sulting from self-avoidance, which is not important for 
strong confinement and in the stiff limit, and is negligi- 
ble even in the flexible limit if the number of segments is 
below N ss 500. At least 10 6 MC steps per segment were 
performed, to obtain averages and statistical errors. Our 
MC procedure was validated by evaluating known quan- 
tities for polymers in bulk [l6| . 

In a typical experimental setup measuring the shape 
of a biopolymer in confined geometry, the filaments are 
labeled with some fluorescent dye and recorded over the 
time resolution window of the camera. This results in an 
intensity profile for the emitted light that corresponds, 
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FIG. 3: Reduced mean-square end-to-end distance A J? 2 as a 
function of the collision parameter c for a series of flexibilities 
e indicated in the graph. The universal scaling curve (solid 
line) asymptotes the analytical result (long-dashed line) in the 
limit of strong confinement. Short-dashed lines are guides to 
the eye for e > 1. Inset: orientational order parameter S as 
a function of c for e = 0.01, 0.1, 1, 10. 



in our theoretical model, to the marginal probability dis- 
tribution function of the positions of the constituent seg- 
ments in a plane containing the tube axis. Figure [5] shows 
this function as obtained from our MC simulations for 
intermediate values of the collision and stiffness parame- 
ters: c = 4 and e = 2. This picture nicely illustrates the 
shape of a biopolymer in confined geometry. Given suffi- 
cient experimental resolution it may even be possible to 
resolve the bimodality of the distribution clearly visible 
in the MC data. 

Beyond such a qualitative impression of the shape 
of the polymer, other more quantitative measures may 
characterize better its conformations. We start our dis- 
cussion with the tangent-tangent correlation function 
4>(s) = (t(s) -t(0.5)), measured from the center of the 
filament. Figure [2] shows the results of our MC simu- 
lations in the stiff regime for s = 0.1 and a range of 
collision parameters c. For small c the data show the 
expected exponential decay of a free filament. For c > 1 
confinement effects become visible and deflections start 
to affect the correlations for distances comparable to the 
Odijk length scale 1/c. For strong confinement, </>(s), af- 
ter an initial exponential decay and a shallow minimum 
at s — 0.5 ~ 1/c, reaches a broad plateau before correla- 
tions decay again in a small boundary layer of size 1/c. 
All these features, but the boundary layer effect, are well 
captured by a formula (solid line in Fig. [2]) easily obtained 
in the WBR approximation from Eq. (H|) : 



(t(a)-t( S ')) = l 
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(5) 



The quantity which best characterizes the elongation 
of the polymer is the mean-square end-to-end distance 
(i? 2 ). As already noted, this is exactly known for uncon- 
fined chains, (-R 2 ) = 2i 2 e~ 2 (e — 1 + e~ e ), but can be 



calculated for strong confinement in the WBR approxi- 
mation. We find: 



L 2 



1 



\/2 e c sin 



(6) 

The validity of this formula, obtained in the WBR limit, 
extends to the case of weak confinement: this result cor- 
rectly recovers the exact result for unconfincd chains in 
the stiff limit. Upon defining the reduced end-to-end dis- 
tance 

AR*(c) = ((R\-(R\)/(L 2 -(r\), (7) 

one finds that this quantity is a function of the collision 
parameter c only. This suggests to look for a data col- 
lapse in the MC data. In fact, as can be inferred from 
Fig. [3J the reduced end-to-end distance is a function of 
the collision parameter only once the flexibility param- 
eter falls below e ~ 1 — i.e., in the stiff regime. This 
implies that there is a single master curve characterizing 
the shape of a stiff polymer. 

The analytical results capture the MC results only in 
the limit of very strong confinement. This is due to 
the fact that the WBR approximation assumes the fil- 
ament to be perfectly aligned with the tube axis, which 
is strictly valid only if c ^> 1 . In the regime of weak con- 
finement the primary effect of the geometric constraints 
is to align the filament with the tube axis. We expect 
this alignment to start once the filament length (more 
precisely the end-to-end distance) becomes comparable 
with the tube diameter — i.e., for c 3 ~ e. To render this 
statement quantitative we define the orientational order 
parameter 
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(8) 



where $ is the angle of R with respect to the tube axis. In 
fact, as can be inferred from the inset of Fig. [31 there is an 
intermediate confinement regime where the onset of ori- 
entational order precedes filament elongation. This has 
important implications for the confinement of biopoly- 
mers like F-actin and microtubules in artificial channels 
and cellular systems. For instance, one estimates that 
for F-Actin with L = 2/im this window in tube dimen- 
sions ranges from d = 2/j.m down to d w 0.4/im. In 
this window the free energy cost for confinement will not 
be given by the Odijk estimate F ~ kgTc but by the 
constraint on the orientational degrees of freedom — i.e., 
F ~ ksTln^L/d). This intermediate regime becomes 
less pronounced with increasing flexibility. Actually, even 
very long (self-avoiding) polymers are known to have an 
instantaneous prolate shape [l7l . [l8| . This anisotropy in 
the radius of gyration tensor is rather due to entropy ef- 
fects [l9| than the energy of bending as for stiff biopoly- 
mers discussed here. Despite the different physical origin 
its consequences are that also flexible polymers in con- 
finement orient first before changing their shape (2(j| . 



4 



1 

0.8 
0.6 
0.4 
0.2 





• 




= 0.1 


X 




= 0.5 


□ 


£ 


= 1 


■ 


£ 


= 5 


G 


£ 


= 10 



0.1 



10 

c/e 



100 



1000 



FIG. 4: Scaling plot for the mean-square end-to-end distance 
(-R 2 ) versus c/e = £ P /Ld for various flexibilities e. Universal 
scaling curve as estimated from the simulation (solid line) and 
obtained in the WBR limit (dashed line). 



For flexible polymers, with e > 1, the initial effect of 
confinement is to increase R\\ but decrease R. There is a 
clear dip in AR 2 , which can be explained by the following 
geometric argument. Consider an initially randomly ori- 
ented end-to-end vector R. Then weak confinement will 
predominantly reduce the magnitude of the component 
of this vector perpendicular to the tube axis but leave the 
parallel component unchanged. This obviously leads to 
a decrease in the magnitude of the end-to-end distance. 
We have checked this argument by measuring the pro- 
jection of the mean-square end-to-end distance onto the 
tube axis (data not shown), which indeed does not show 
an initial decrease but increases monotonically. A dip in 
the radius of gyration Rq (not the end-to-end distance) 
has previously been reported for very long (e> 1) self- 
avoiding polymers [2(| ■ There it was argued to be a con- 
sequence of an intermediate regime, where the principal 
components of the radius of gyration tensor are reduced 
before confinement excluded- volume effects lead to an in- 
crease in Re- 
in the limit of strong confinement, where the Odijk 
deflection length is either much smaller than the total 
length (for stiff polymers) or much smaller than the per- 
sistence length (for flexible polymers), one expects (R 2 ) 
to become independent of the filament length. In order 
to show this explicitly we have replotted our MC data in 



Fig. [4] as a function of the number of collisions within a 
persistence length, c/e = l p /L^. Indeed, all curves for 
different flexibilities merge into a master curve for strong 
confinement. In contrast to the master curve for weak 
confinement, which applies only for stiff chains, this mas- 
ter curve captures all chain flexibilities. The numerical 
result asymptotes the analytical formula obtained in the 
WBR limit, but only for c/e quite large. This shows that 
for intermediate confinement the orientation and the less 
constrained ends of the filament contribute significantly 
to the conformations; both of these effects are neglected 
in the WBR limit. 

Advances in microfabrication and nanofabrication 
technologies have made it possible to confine biopolymers 
to topographical structures whose dimensions arc com- 
parable to or even smaller than their persistence length. 
This opens a range of novel possibilities to visualize and 
manipulate DNA and cytoskeletal filaments. Here, we 
have presented an extensive theoretical analysis of the 
shape and conformations of biopolymers resulting from 
strong confinement and identified and quantified a range 
of novel scaling regimes. We make specific predictions for 
experimentally accessible quantities like the density pro- 
file or the orientation and apparent length of a biopoly- 
mer in a channel. Our calculations may provide a road 
map for a clear identification of the possible scaling sce- 
narios involved in the manufacture of nanofluidic and 
microfluidic devices. At the same time, our analysis is 
a first step towards a quantitative master curve connect- 
ing the apparent length and the actual length of DNA 
in nanochanncls, which has important implications for 
experimental realizations aimed at a rapid screening of 
entire genomes in contrast to gel electrophoresis. Finally, 
they shed some light on the effect of cellular crowding on 
the conformation of cytoskeletal filaments. 
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